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Abstract 

It is proved that feedback classification of a linear system over a commutative 
von Neumann regular ring R can be reduced to the classification of a finite family 
of systems, each of which is properly split into a reachable and a non-reachable 
part, where the reachable part is in a Brunovski-type canonical form, while the 
non-reachable part can only be altered by similarity. If a canonical form is known 
for similarity of matrices over R, then it can be used to construct a canonical form 
for feedback equivalence. An explicit algorithm is given to obtain the canonical 
form in a computable context together with an example over a finite ring. 



1 Introduction and motivation 

Let R be a commutative ring with 1. An m-input, n- dimensional linear system over R is a 
pair of matrices (A, B), with A e R nxn , B e R nxm . The control theoretic background for 
this terminology is nicely described in [1] Introduction] . Within this context two systems 

(A, B) and (A', B') are called feedback equivalent (f.e., rJ) if there exist matrices P, Q, K 
of suitable sizes such that (A', B') = (P(A + BK)P~ l , PBQ), with P, Q invertible. 

Note that equivalence and similarity of matrices over R are particular cases of feed- 
back equivalence, indeed, if (A, B) and (A', B') are feedback equivalent, then B, B' are 
equivalent, while the equivalence of (A, 0), (A',0) implies similarity of A, A'. 

Classification and canonical forms under feedback equivalence for linear systems are 
classical problems in linear control theory. While a general solution for arbitrary rings 
can not be expected several specific results have been obtained for special classes of 
commutative rings. If R is a field, controllability indices and Brunovski canonical form 
for reachable systems together with Kalman's decomposition give a complete description. 

For larger classes of rings a variety of partial results have been obtained so far, e.g. a 
partial canonical form for reachable systems over principal ideal domains [5] or a canonical 
form for weakly reachable single-input systems over Bezout domains [5]. 

In |3], [H] it was observed, that certain basic control problems have a solution if and 
only if the underlying ring is regular (in the sense of von Neumann, see definition below): 

- local and global feedback equivalence over R are equivalent iff R is regular. 

- cyclizability of the reachable part of systems is always possible iff R is regular. 
For details see [U Theorem 4.1] and [91 Theorem 4]. 
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At the same time regular rings have nice linear algebraic properties, for example: 

- R nxn is regular iff R is regular [5J. 

- The tensor criterion for matrix similarity of Byrnes and Gauger is valid iff the underlying 
ring is regular [7]. 

This encourages to study feedback classification and canonical forms over the large class 
of commutative regular rings. A commutative ring R is said to be (von Neumann) regular 
if for any a in R there exists i6i? such that a 2 x = a. 

The following mostly well-known properties of regular rings will be of use for us later, 
see [6J for details. 

- Any finitely generated ideal is principal (i.e. R is a Bezout ring) and generated by 
an idempotent. 

- Any element is the product of a unit with an idempotent. 

- R has stable range 1, which for a Bezout ring means that whenever aR + bR = dR 
for some a,b,d G R then there exists c G R such that a + be = ud, with u a unit. 

- R is an elementary divisor ring, which means that there exists a 'Smith Form': every 
matrix B is equivalent to a diagonal matrix with diagonal entries d\\d2\- • - \d r . In 
fact, the Smith Form is unique, if the d^s are idempotents (see Lemma H] below) . 

- Any finitely generated sub module of R n has a finitely generated complement. 

Note also that any field is a regular ring and more generally any subring of any product 
of fields and in particular many rings of functions comprising many finite rings. 

2 Preliminaries 

For an n-dimensional system E = (A, B) over a commutative ring the right image of 
[B, AB, . . . , A n ~ 1 B] is called the module of reachable states. £ is reachable if TZs = R n - 
When R is regular, then TZy, is always a (finitely genarated) direct summand of R n with 
a finitely generated complement. If TZy, has a finite basis which can be completed to a 
basis of R n , then the system £ is f.e. to a system of the form 



'A' 


A"' 




~B'~ 





A'" 
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where (A', B') is reachable. This is usually called a Kalman decomposition, which we shall 
call strong in case A" = 0. The reachable system (A, B) is in so called Brunovski form if 

A = diag(A, ...,A,),B= [diagOBi, ...,B 8 ),0] 

where Ai = [e^ 2 \ . . . , e^ ni \ 0] , B{ — &^ with the j-th canonical basis vector in R ni for 
1 < i < s and where ri\ > • • • > n s , n\ + • • ■ + n s = n. The rij are called controllability 
or also Kronecker indices. They are a complete set of invariants under feedback in case 
R is a field. See e.g. [TUl chapters 3.3 and 4.2 ] 

Typically regular rings contain many idempotent elements. Therefore we now collect 
some of their basic properties. 

Lemma 1 Let R be a commutative ring and e\, . . . , pairwise orthogonal idempotents 
(titj = ifi^j) and such that Yli=i e i — 1- Then: 

R = e x R © • • • © e k R 

Note that the ideal CiR is at the same time a regular ring with as unit element. 



2 



An immediate consequence of such a decomposition is, that any matricial equation 
M = M' holds over R if and only if ejM = e^M' holds for all i, in particular: 

Lemma 2 -A matrix P is invertible over R iff eiP is invertible over eiR for all i. 

- Two matrices A, A' are similar (resp. equivalent) over R iff CiA^CiA' are similar 
(resp. equivalent) over CiR for all i. 

- Two systems {A,B) and (A',B') are f.e. over R iff (eiA,eiB) and (eiA',eiB') are 
f.e. over eiR for all i. 



Next, we recall basic facts about feedback equivalence. 



Lemma 3 Let R be a commutative ring with 1, and consider the system (A, B) over R 
of size (n, m), < m < n given by 



A 




Bx A, 



B 



with I G R mxm an % 
appropriate sizes. Then, one has: 



identity matrix, A\ e R^ n m ) x ( n m ) and the remaining blocks of 



(i) The feedback equivalence class of(A,B) is uniquely determined by the feedback equiv- 
alence class of the (n — m,m) system (Ai,B{) in the following sense: two systems 



( 




B x A x 



Y ( 




B[ A[ 



) 



are equivalent if and only if (A\, Bi), (A[, B[) are equivalent, 
(ii) If R is a principal ideal domain or has stable range 1 (in particular, if R is regular), 
(i) is also valid if B has any number of additional zero columns, i.e. it has the block 
form: 

\I Ol 



B 







(Hi) If (Ai,Bx) is a reachable system in Brunovski canonical form then also (A,B) 
is reachable and can be transformed just by permutation of rows and columns, to 
Brunovski form. Moreover the controllability indices Ki > ■ • • > n m , Ai > • • • > A m 
(possibly some of the Ki have to be put to zero) of (Ai,Bi) and (A,B) are related 
by Aj = ^ + 1 for 1 < i < m. 



Proof, (i) The proof is straightforward and can be adapted from [3 Lemmas 2.1 and 
2.2]. 

(ii) See [SJ Propositions 2.4 and 2.5]. 

(iii) When (A\, B\) is in Brunovski form, then for some s we have 

Ae {1) = e (m+1) , A 2 e {1) = e (m+2) , . . . , A s e {1) = e {m+s) . 

Similar relations are valid for e^ 2 \ . . . , e^ m \ Thus an appropriate renumbering of the 
canonical basis vectors gives a Brunovski form and the relation for the indices is straight- 
forward. 

□ 
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3 Feedback classification and canonical forms 



Let R be a regular ring, and consider a system E = (A, B) of size (n, m) over R. We 
will always tacitely assume B ^ 0. In the sequel we will construct recursively a normal 
form and a complete set of invariants for the feedback classification of E. The following 
reduction step will be essential: 



Lemma 4 By operations of the feedback group any system (A, B) can be assumed to have 
the following form 

A=(a ij ) n i j=v B = 



with d\\d s \ ■ ■ ■ \d r , where all the di 's are nonzero idempotents and where for i = 1, . . . , r, 
the i-th row of A is orthogonal to di. Furthermore, the elements d\, . . . ,d r are invariant 
under feedback equivalence. If r = n or r = m the corresponding 0-blocks will not occur. 



di 









d r 











Proof. Since R is an elementary divisor ring and any element of R is a product of an 
idempotent with a unit, the form of B is obtained via matrix equivalence. Moreover, if 
we denote by the element in position (i,j) of A, with a suitable feedback action it can 
be replaced by — diaij = (1 — d^a^, which is orthogonal to di. 

For each i = 1, . . . , r, the ideal of R generated by all the i x i minors of B is invariant 
under equivalence and thus invariant under feedback equivalence of the system (A, B). 
But this ideal is clearly generated by d\ ■ ■ ■ di — di, and an ideal in R cannot be generated 
by two distinct idempotents. □ 

We are now ready to solve the feedback classification problem for systems over a von 
Neumann regular ring R. A first attempt in this direction could be to reduce the problem 
to the classification of systems over all residue fields (see [H Theorem 4.1]): two systems 
over R are equivalent if and only if they are equivalent over the residue field R/xn, for 
all maximal ideals m. However, there may exist an infinite number of maximal ideals in 
the regular ring R. But, as we will prove in the next theorem, the feedback classification 
problem of a system over R can be reduced to the classification of a finite family of systems 
which behave like systems over fields. 



Theorem 5 (Canonical forms and invariants) Let E = (A, B) be a system of size 
(n,m) over a regular ring R. Then: 

(i) There exists a finite family of idempotents {e^}, pairwise orthogonal and with sum 
1, such that for each k, the system (e k A, e k B) over the regular ring e k R is feedback 
equivalent to a system in strong Kalman decomposition, with the block form 



A k 




Bk 


C k _ 








)■ 



where the pair (A k ,Bk) is reachable and in Brunovski canonical form over e k R. 
(ii) With these notations, a complete set of feedback invariants for (A,B) consists of: 
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- The idempotents {e^}, which can be obtained successively from the invariant 
factors of a finite number of matrix equivalences. 

- The similarity classes of the matrices C k . 

- The controllability indices of the systems (A k ,B k ). 



Proof. We will proceed by induction on n. 

• Proof of (i) and (ii) in the case n — 1. We can assume without loss of generality that 
the system E is given in the form ([1]): A = [a],B = [d, 0-- -0], with d idempotent and 
ad = 0. Let us consider the partition R = dR © (1 — d)R, and build the systems: 

(l-tf)£= ([(1 -d)a],{0 ■■■()]) and rfS = ([0], [d, • • ■ 0]) = d([0], [1, 0, . . . , 0]) 

It is clear that both systems are decomposed in strong Kalman form, the first one without 
reachable part, and the second one already in Brunovski form over the ring dR, in which 
d is the 1-element. Moreover, the element d is obtained from the equivalence of B and 
is invariant by Lemma HJ making 1 — d also invariant. Now, by Lemma El the feedback 
equivalence class of £ is completely determined by that of <i£ (with a single invariant 
d), and that of (1 — where the first entry (1 — d)a in principle can be reduced by 
similarity, but in the 1 x 1-case this has no effect. 

• Proof of (i) for n > 1. Let E be given as in ([1]), and consider the following family of 
idempotents (pairwise orthogonal and with sum 1): 

e = 1 — d\ 
e± = di(l - d 2 ) 

e r _i = d r _i(l — d r ) 



For % = 0, we have that eo = 1 — d\ is orthogonal to d±, and hence also orthogonal to all 
the remaining d^s, so that e$B = 0, which means that eoE trivially is a strong Kalman 
decomposition with no reachable part. This step is omitted if di — 1. 

For i = 1, . . . , min{r, n — 1}, note that is a multiple of di, . . . , di and orthogonal to 
di + i, . . . , d r and to the rows 1, . . . , i of A. Note also, that = if di = dt+i- 

Therefore we obtain for 1 < i < min{r, n — 1} and if e-i ^ 0: 



e^S = (eiA, eiB) = ( 











ej 


" 


e-iBi 


&iAi 











) 



with Ai G R^ n *- )x( ' n l \ I an i x i identity block, and the remaining blocks of appropriate 
sizes. 

Since (ejAj, ej-Bj) is of size (n — i, i) (with n — i < n), by the induction assumption one 
can assume that there exists a finite partition e^i? = ® . e^R, such that for each j, the 
system (e^Ai, e^Bi) over the ring e^R is equivalent to one of the form 



) 



Aij 




Bij 


C i:i _ 








with (Aij,Bij) reachable and in Brunovski canonical form. 



5 



In virtue of Lemma[3],(i),(ii), we obtain 














/ e / 

CijH = [SijA, eijB) ~ ( 




















Cijl 


" 















The system on the right side can also be written as 

( 












Bij 













Cij _ 



tijl 

_0 0_ 





By Lemma [31 (iii) we know that 



" " 




djl 


R A-- 
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with (Aij, By) in Brunovski form. If we let = C^, then we obtain finally 



A^ 




Bij 


_ dj_ 




_ 



with the desired properties for the righthand side system. This completes the study of 
the systems e^E, for % — 1, . . . , min{r, n — 1}. 

If r = n, the system e n S is already in Brunovski form, because e n A is zero and e n B = 
[e n /|0], with I an n x n identity block. 

To sum up, we have obtained a finite partition of R with idempotents {e^} such that 
for each the system e^Y, is a strong Kalman decomposition with reachable part in 
Brunovski form, as required. 

• Proof of (ii) for n > 1. With the preceeding notations, the elements {e^} are con- 
structed from the Smith form of B and satisfy the conditions of Lemma [TJ from which 
it follows by Lemma [2] that the classification of £ over R is reduced to the classification 
of the systems e^S over the regular rings CiR (at this step we need only consider those 
indices % for which 7^ 0). A complete set of invariants for S will be given by the elements 
3q, . . . , e r , together with all the invariants obtained recursively. 

For % = we have seen in (i) that e^B = 0, i.e. the feedback class of (egA, cqB) is 
reduced to the similarity class of e$A. 

For 2 = 1,..., min{r, n — 1}, by Lemma [3]the feedback class of e^S is completely deter- 
mined by the feedback class of the system {eiAi, tiBi), which by the induction hypothesis 
(see the notation in the proof of part (i)) reduces to a collection of idempotents {e^} 
obtained from equivalences of matrices, the similarity class of all matrices C^, and the 
controllability indices of all the systems (Aij,Bij), which are directly related to the con- 
trollability indices of the augmented systems (A^, B^) as was explaineed in Lemma [31 (iii) • 

Finally, if r = n, the additional system e n S has the trivial form (0, [e n J|0]) and does 
not provide any new invariant. This completes the proof. □ 
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4 The reachable and the single input case 



As an immediate consequence, reachable systems can be completely described without 
reference to similarity, i.e. only with idempotents and Brunovski blocks. 

Corollary 6 (Reachable systems) With the above notations, every reachable system 
over a regular ring is equivalent to a direct sum of systems in Brunovski canonical 
form, and the feedback equivalence class can be described by a family of idempotents and 
Brunovski indices. 

Proof. With the notations of the previous theorem, reachability of a given system (A, B) 
is equivalent to reachability of all the systems 



(e k A, e k B) ~ ( 



A k 






C k _ 








)- 



which is only possible if none of the blocks C k appear, i.e. one can remove all references 
to similarity in the previous theorem. □ 

The description via a finite collection of idempotents seems essential and unavoidable 
for an explicit exposure of some kind of normal form. Nevertheless, as in the classical 
field case, over a regular ring it is possible to classify reachable linear systems in a more 
global, but at the same time more abstract way, by the following sequence of submodules 
of the module TZy, of reachable states. 

For an n-dimensional system E = (A, B) over a regular ring R let 

N% =im[B\AB\---\A k - 1 B], 

for k — 1, . . . , n. 

Proposition 7 (Reachable systems and the modules Let E = {A,B) and 
E' = (A 1 , B') be two reachable systems of size (n,m) over a regular ring R. Then, the 
following statements are equivalent: 

(i) E and E' are feedback equivalent. 

(ii) The R n -submodules and Nj? are isomorphic for all r — 1, . . . , n. 

Proof, (i) =>■ (ii) is trivial, (ii) =>- (i) can be derived directly from results in [I]. Let m 
be an arbitrary maximal ideal of R. If we denote by S(m) and S'(m) the extensions of 
E, E' to the residue field R/ra (by reduction modulo m), since = we must have 

k k 

for k — 1, . . . , n. But the dimensions of the R/m- vector spaces N^ m \ ^ characterize 
completely the equivalence classes of the reachable systems E(m),E'(m) over the field 

R/m (see [U Proposition 2.5]), therefore it follows that S(m) ~ S'(m) over R/m. Now, 

by [H Theorem 4.1] we can conclude that E ~ E' over the regular ring R. □ 

Proposition [7] is yet another instance where regular rings behave "classically". Even 
for principal ideal domains a corresponding general result is not possible because of a 
counterexample in [21 section 6.2]. 

As a consequence of Proposition [7] we obtain 
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Corollary 8 The feedback equivalence class of a reachable system £ = (A, B) of size 
(n, m) over a regular ring R is completely determined by the invariant factors of the 
matrices [B\AB\ ■ ■ ■ \A k ~ 1 B], for k — 1, . . . , n. 



Proof. The result is clear, since for each k = l,...,n, the columns of the matrix 
[£?|AB| • • • generate the R- module N^, whose isomorphism class is determined 

by the invariant factors of the generating matrix. □ 

In the single-input case, the reachable part of a system can be transformed to a gener- 
alized controller canonical form. All notations before. 



Proposition 9 (Single input systems) . 

(i) A canonical form for a system £ = (A, B), B ^ 0, of size (n, 1) over a regular ring 
R is given by 



A 






•• 


d 2 


o '• 





d 3 ' ■ 





... o 











A* , B 



'd 1 








with idempotents di\d 2 \ ■ ■ ■ \d n , possibly d r+ i 
with 

J. ~ TO 

A = A 1 . dl 



d n = for some r > 1, and 



A dl( i_ d2 ) 



+ 





A 





<i r — i(l— d r ). 



where A di(1 _ di+l) e J R(™- i ) x ("- i ) for 1 < i < r and d i+1 A di{l „ di+l) = forO<i<r. 
(ii) For k = 1, . . . , n one has 



n (a,b) ^ n (a,b) = die (i) R + ... + d k e {k) R = d 1 R®d 2 R®---® d k R 



(Hi) The elements d±, . . . , d n can be obtained as the idempotent invariant factors appear- 
ing in the Smith form of the matrix [B, AB, . . . , A n ~ 1 B]. 

Proof. 1.: For a system E = (A,B) of size (n, 1), the reduced form ([T]) yields only 
one idempotent invariant factor d\, i.e. two idempotents eo = 1 — d\ and e\ = di as in 
the proof of Theorem [51 The system eoS can only be reduced further by similarity and 
is characterized by some corresponding canonical form Ai_ dl , while the class of d{L is 
reduced to that of a system (d\Ai, diBi) of size (n—1, 1). Continuing this iteration r times, 
at each step i we will have a system {diA^ diBj) of size (n — i, 1), with <ij +1 the idempotent 
invariant factor of diBi (di\d i+ i), and A d .^ di+l ) a canonical form of — d i+ i)Ai for 
similarity over the ring di(l — di + \)R. Multiple use of Lemma El(i) leads to the stated 
result. 

2. and 3.: The modules Nj? are invariant under feedback and remain isomorphic under 
feedback equivalence. Therefore 2. and 3. are straightforward consequences of 1. □ 
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5 Implementation 



The procedure given in Theorem[5]is constructive, which gives rise to an explicit algorithm, 
provided Smith normal forms are computable over the ring R, and a canonical form is 
known for the similarity relation. 

Algorithm 10 (Canonical forms) . 

INPUT: (e,A,B), whith (A,B) a system and e an idempotent (initially e = 1) 

OUTPUT: A list of lists [e^, Ak, Bk], where ek is an idempotent and (Ak, Bk) is a canonical 
form of (A, B) over the ring CkR- 

1. Replace (A,B) by its reduced form (TJP, with d\\d2\ ■ ■ ■ \d r idempotents. 

2. Initialize L=[ ] (this will be the output). 

3. Define {eo = e — d±, e% = di(e — d-i), . . . , e r _i = d r -i(e — d r ), e r = d r }. 
4- For each ^ 0, do steps 5-8. 

5. If e-iB = 0, add [e», Aj, 0] to L, where Ai is a canonical form of A for similarity over 
CiR. Proceed with next ej. 

6. If i = r = n, add [e^, 0, CiE\ to L. 





7. Extract blocks eiA iy eiBi from CiA 



&%T$i ^iAi 



8. Recursive call with input (ej, e%Ai, eiBj) and output a list 



&ij j Aij , Bij 



J .I 



9. For each j , add to L : 



R A ■■ 



Bij I 





10. If ei — for all i (L is still [ ]), output [[e, A, 5]] , otherwise output L. 

Example 11 We have used a symbolic software to compute canonical forms for rings 
Z/(oK), where d is a squarefree integer. The following system over the ring Z/(210Z) 
splits into 3 subsystems: 



105S 



S 



105 105 105 
105 105 
105 
105 105 105 



36S 



21 158 169 147 
138 208 43 135 
67 46 190 100 
167 36 81 203 







36 



70S 



178 152 60 58 ' 
90 186 36 120 
102 96 30 198 
140 40 42 146 

70 
70 
140 
36 0" 
36 
36 




70 






We see that the classification of (A, B) reduces to the similarity classification of a 4 x 4 
block of the first system, and the lxl block 140 appearing in the uncontrollable part 
of the second system, while the third system is reachable and completely described by 
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the sizes of the idempotent and zero blocks which determine its canonical form. E.g, the 
Brunovski form of 36S is 



36 ■ ( 



"0 





0" 




"1 








0" 


1 




































1 
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6 Concluding Remarks 

Our results show that over a regular ring feedback classification of systems still behaves 
quite similar to the classical context over a field. But typically one has to handle a 
finite collection of systems at the same time. Each of these "parallel" systems can be 
transformed constructively into a strong Kalman decomposition, where the reachable 
part is in Brunovski canonical form and where the non reachable part can only be altered 
further via matrix similarity. Only if a normal form for similarity of matrices is known, 
as is the case over fields, our results give a complete normal form. 

The connection of our normal form to the elementary divisors of the reachability sub- 
module, which we state in Proposition [9] has not yet been worked out in the multi-input 
case. 
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